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We study the non-abelian statistics characterizing systems where counter-propagating gapless 
modes on the edges of fractional quantum Hall states are gapped by proximity-coupling to su- 
perconductors and ferromagnets. The most transparent example is that of a fractional quantum 
spin Hall state, in which electrons of one spin direction occupy a fractional quantum Hall state of 
v = 1/m, while electrons of the opposite spin occupy a similar state with v = —1/m. However, we 
also propose other examples of such systems, which are easier to realize experimentally. We find 
that each interface between a region on the edge coupled to a superconductor and a region coupled 
to a ferromagnet corresponds to a non-abelian anyon of quantum dimension y/2m. We calculate the 
unitary transformations that are associated with braiding of these anyons, and show that they are 
able to realize a richer set of non-abelian representations of the braid group than the set realized 
by non-abelian anyons based on Majorana fermions. We carry out this calculation both explicitly 
and by applying general considerations. Finally, we show that topological manipulations with these 
anyons cannot realize universal quantum computation. 



I. INTRODUCTION 

Recent years have witnessed an extensive search for 
electronic systems in which excitations ("quasi-particles") 
follow non-abelian quantum statistics. In such sys- 
tems, the presence of quasi-particles, also known as "non- 
abelian anyons'S^, makes the ground state degener- 
ate. A mutual adiabatic interchange of quasi-particles' 
positions 4 implements a unitary transformation that op- 
erates within the subspace of ground states, and shifts 
the system from one ground state to another. Remark- 
ably, this unitary transformation depends only on the 
topology of the interchange, and is insensitive to im- 
precision and noise. These properties make non-abelian 
anyons a test-ground for the idea of topological quan- 
tum computation 5 . The search for non-abelian systems 
originated from the Moore-Read theory^ for the v = 5/2 
Fractional Quantum Hall (FQH) state, and went on to 
consider other quantum Hall states^! spin systems^, 
p-wave superconductors^^ 2 !, topologi cal in sulators in 
proximity coupling to superconductors^^ and hybrid 
systems of superconductors coupled to semiconductors 
where spin-orbit coupling is strongEaE] Signatures of 
Majorana zero modes may have been observed in recent 
experiments^ 2 ** 2 ^. 

In the realizations based on superconductors, whether 
directly or by proximity, the non-abelian statistics results 
from the occurrence of zero-energy Majorana fermions 



bound to the core s of vor tices or to the ends of one di- 
mensional wires^** 21 * 27 * 28 !. Majorana-based non-abelian 
statistics is, from the theory side, the most solid pre- 
diction for the occurrence of non-abelian statistics, since 
it is primarily based on the well tested BCS mean field 
theory of superconductivity. Moreover, from the experi- 
mental side it is the easiest realization to observe 22 !. The 
set of unitary transformations that may be carried out 
on Majorana-based systems is, however, rather limited, 
and does n ot all ow for universal topological quantum 
computation^^]. 

In this work we introduce and analyze a non-abelian 
system that is based on proximity coupling to a su- 
perconductor but goes beyond the Majorana fermion 
paradigm. The system we analyze is based on the 
proximity-coupling of fractional quantum Hall systems 
or fractional quantum spin Hall systems^ to supercon- 
ductors and ferromagnetic insulators (we will use the 
terms "fractional topological insulators" and "fractional 
quantum spin Hall states" interchangeably). The start- 
ing point of our approach is the following observation, 
made by Fu and Kane^when considering the edge states 
of two-dimensional (2D) topological insulators of non- 
interacting electrons, of which the integer quantum spin 
Hall stat d 32 * 33 -* is a particular example: In a 2D topo- 
logical insulator, the gapless edge modes may be gapped 
either by breaking time reversal symmetry or by breaking 
charge conservation along the edge. The former may be 
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broken by proximity coupling to a ferromagnet . while the 
latter may be broken by proximity coupling to a super- 
conductor. Remarkably, there must be a single Majorana 
mode localized at each interface between a region where 
the edge modes are gapped by a superconductor to a re- 
gion where the edge modes are gapped by a ferromagnet. 

Our focus is on similar situations in cases where the 
gapless edge modes are of fractional nature. We find 
that under these circumstances, the Majorana operators 
carried by the interfaces in the integer case are replaced 
by "fractional Majorana operators" whose properties we 
study. 

We consider three types of physical systems. The first 
(shown schematically in Fig. [TJl) is that of a 2D frac- 
tional topological insulatorSS, that may be viewed as a 
2D system in which electrons of spin-up form an FQH 
state of a LaughlirPSl fraction v = 1/m , with m being 
an odd integer, and electrons of spin-down form an FQH 
state of a Laughlin fraction v = —1/m. 

The second system (|T|d) is a Laughlin FQH droplet 
of v = 1/m, divided by a thin insulating barrier into 
an inner disk and an outer annulus. On the inner disk, 
the electronic spins are polarized parallel to the magnetic 
field (spin-up), and on the annulus the electronic spins 
are polarized anti-parallel to the magnetic field (spin- 
down). Consequently, two edge modes flow on the two 
sides of the barrier, with opposite spins and opposite ve- 
locities. Such a state may be created under circumstances 
where the sign of the g-factor is made to vary across the 
barrier. 

The third system is an electron-hole bi-layer subjected 
to a perpendicular magnetic field, in which one layer is 
tuned to an electron spin-polarized filling factor of v = 
1/m, and the other to a hole spin-polarized v = 1/m 
state. In particular, this may be realized in a material 




Figure 1. Schematic setup, (a) A fractional topological 
insulator (FTI) realization. A FTI droplet with an odd filling 
factor 1/m is proximity coupled to ferromagnets (FM) and 
to superconductors (SC), which gap out its edge modes. The 
interfaces between the SC and FM segments on the edge of the 
FTI are marked by red stars, (b) A fractional quantum Hall 
(FQH) realization. A FQH droplet with filling factor 1/m is 
separated by a thin barrier into two pieces: an inner disk, 
and an outer annulus. On either side of the barrier, there are 
counter-propagating edge states, which are proximity coupled 
to superconductors and ferromagnets. 



with a spectrum that is electron-hole symmetric, such as 
graphene. 

In all these cases, the gapless edge mode may be 
gapped either by proximity coupling to a superconductor 
or by proximity coupling to a ferromagnet. We imagine 
that the edge region is divided into 2N segments, where 
the superconducting segments are all proximity coupled 
to the same bulk superconductor, and the ferromagnetic 
segments are all proximity coupled to the same ferromeg- 
nct. The length of each segment is large compared to the 
microscopic lengths, so that tunneling between neighbor- 
ing SC-FM interfaces is suppressed. We consider the 
proximity interactions of the segments with the super- 
conductor and the ferromegnet to be strong. 

The questions we ask ourselves are motivated by the 
analogy with the non-interacting systems of Majorana 
fermions: what is the degeneracy of the ground state? 
Is this degeneracy topologically protected? What is the 
nature of the degenerate ground states? And how can one 
manipulate the system such that it evolves, in a protected 
way, between different ground states? 

The structure of the paper is as follows: In Sec. [IT] we 
give the physical picture that we developed, and summa- 
rize our results. In Sec. IIIII we define the Hamiltonian 
of the system. In S ec. |IV| we calculate the ground state 
degeneracy. In Sec. |V| we define the operators that are 
localized at the interfaces, and act on the zero energy 
subspace. In Sec. |VI| we calculate in detail the unitary 
transformation that corresponds to a braid operation. In 
Sec. |VII| we show how this transformation may be de- 
duced from general considerations, bypassing the need for 
a detailed calculation. In Sec. IVIIll we discuss several as- 
pects of the fractionalized Majorana operators, and their 
suitability for topological quantum computation. Sec. |IX| 
contains some concluding remarks. The paper is followed 
by appendices which discuss several technical details. 



II. THE PHYSICAL PICTURE AND SUMMARY 
OF THE RESULTS 

There are three types of regions in the systems we con- 
sider: the bulk, the parts of the edge that are proximity- 
coupled to a superconductor, and the parts of the edge 
that are proximity-coupled to a ferromagnet. 

The bulk is either a fractional quantum Hall state or 
a fractional quantum spin Hall state. In both cases it 
is gapped and incompressible, and its elementary excita- 
tions are localized quasi-particles whose charge is a mul- 
tiple of e* — e/m electron charges. In our analysis we 
will assume that the area enclosed by the edge modes en- 
closes quasi-particles of spin-up and quasi-particles 
of spin down. These quasi-particles are assumed to be 
immobile. 

In the parts of the edge that are coupled to a super- 
conductor the charge is defined only modulo 2e, because 
Cooper-pairs may be exchanged with the superconduc- 
tor. Thus, the proper operator to describe the charge on 
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a region of this type is e MrC " , with Qi being the charge 
in i'th superconducting region. Since the superconduct- 
ing region may exchange e* charges with the bulk, these 
operators may take the values e l7T <a/ m j with an integer 
whose value is between zero and 2m — 1. The pairing 
interaction leads to a ground state that is a spin singlet, 
and thus the expectation value of the spin within each 
superconducting region vanishes. As we show below, the 
Hamiltonian of the system commutes with the operators 
e "Q; m the limit we consider. For the familiar m — 1 
case, these operators measure the parity of the number 
of electrons within each superconducting region. 

The edge regions that are proximity-coupled to ferro- 
magnets are, in some sense, the dual of the superconduct- 
ing regions. The ferromagnet introduces back-scattering 
between the two counter-propagating edge modes, lead- 
ing to the formation of an energy gap. If the chemical 
potential lies within this gap, the region becomes insulat- 
ing and incompressible. Consequently, the charge in the 
region does not fluctuate, and its value may be defined 
as zero. The spin, on the other hand, does fluctuate. 
Since the back-scattering from spin up electron to spin 
down electron changes the total spin of the region by two 
(where the electronic spin is defined as one unit of spin) , 
the operator that may be expected to have an expecta- 
tion value within the ground state is e l7rSi , with §i being 
the total spin in the i'th ferromagnet region. Again, spins 
of 1/m may be exchanged with the bulk, and thus these 
operators may take the eigenvalues e l7rSi / m , with s an in- 
teger between zero and 2m — 1. The Hamiltonian of the 
system commutes with the operators e * in the limit 
we consider. 

The operators e l ' K< ^ i and e mSi label the different do- 
mains in the system, as indicated in Fig |2^i. They satisfy 
a constraint dictated by the state of the bulk, 
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For the familiar m = 1 case there are only two possible 
solutions for these constraints, corresponding to the two 
right hand sides of Eq. ([lj being both +1 or both — 1. 
For a general m, the number of topologically distinct 
constraints is 2m 2 , since Eqs. are invariant under 
the transformation where n+ —¥ n+ ± m together with 
n\. — > nj. ± m. These sets may be spanned by the values 
< < 2m — 1 and < nj. < m — 1 . 

The degeneracy of the ground state may be understood 
by examining the algebra constructed by the operators 
gijrQi an( j giTrS^ ^ g we snow m tfie next section, the 
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where, in the last equation, 1 < j,l < N (see Fig [2^, for 
the enumeration convention) . As manifested by Eqs. Q, 
the pairs of operators e 1 *^ , e"^ k=1 s *> form N — 1 pairs 
of degrees of freedom, where members of different pairs 
commute with one another. It is the relation between 
members of the same pairs, expressed in Eq. (J2|, from 
which the ground state degeneracy may be easily read 
out. As is evident from this equation, if is a ground 

state of the system which is also an eigenstate of e I7r( ^ , 
then 2m — 1 additional ground are ^e I7r S=i 

where k is an integer between 1 and 2m — 1 . With N — 1 
mutually independent pairs, we reach the conclusion that 
the ground state degeneracy, for a given value of n^,n^, 
is (2m) N - 1 . 

The operators acting within a sector of given rif, 
ni of the ground state subspace are represented by 
(2m) N - 1 x (2m) N - 1 matrices. They may be expressed in 
terms of sums and products of the operators appearing 
in ([2]). The physical operations described by the oper- 
ators e l7r ^ i , e l7r( ^ i can also be read off the relations Q . 
The operator e 47rSi transfers a quasi-particle of charge 
e/m from the i — l'th superconductor to the i'th su- 
perconductor. Since the spin within the superconductor 
vanishes, there is no distinction, within the ground state 
manifold, between the possible spin states of the trans- 
ferred quasi-particle. In contrast, the operator e l * a Q i 
transfers a quasi-particle of spin a = ±1 across the i'th 
superconductor . 

For m — 1 the operators e™ Si and e t7r< ^ i , measuring 
the parity of the spin and the charge in the i'th ferro- 
magnetic and superconducting region, respectively, may 
be expressed in terms of Majorana operators that reside 
at the interfaces bordering that region. A similar repre- 
sentation exists also in the case of m ^ 1. Its details are 
given in Section [Vj 

We stress that the ground state degeneracy is topologi- 
cal, in the sense that no measurement of a local operator 
can determine the state of the system within the ground 
state subspace. For m — 1, this corresponds to the well- 
known "topological protecti on" o f the ground state sub- 
space of Majorana fermionj 29 ! 35 !, as long as single elec- 
tron tunneling is forbidden either between the different 
Majorana modes, or between the Majorana modes and 
the external world. In the fractional case the states in 
the ground state manifold can be labelled by the frac- 
tional part of the spin or charge of the FM/SC segments, 
respectively. These clearly cannot be measured locally. 
Moreover, they can only change by tunneling fractional 
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quasi-particles between different segments; even tunnel- 
ing electrons from the outside environment cannot split 
the degeneracy completely, because it can only change 
the charge and spin of the system by integers. 

Topological manipulations of non-abelian anyons con- 
fined to one dimension are somewhat more complicated 
than those carried out in two dimensions. The simplest 
manipulation does not involve any motion of the anyons, 
but rather involves either a 2tt twist of the order param- 
eter of the superconductor coupled to one or several su- 
perconducting segments, or a 2ir rotation of the direction 
of the magnetization of the ferromagnet coupled to the 
insulating segments 36 . When a vortex encircles the i'th 
superconducting region it leads to the accumulation of a 
Berry phase of 27r multiplied by the number of Cooper- 
pairs it encircles. In the problem we consider, this phase 
amounts to e M( ^ i , and that is the unitary transformation 
applied by such rotation. As explained above, this trans- 
formation transfers a spin of 1/m between the two fer- 
romagnetic regions with which the superconductor bor- 
ders. Similarly, a rotation of the magnetization in the 
ferromagnetic region leads to a transfer of a charge of 
e/m between the two superconductors with which the 
ferromagnet borders. 

A more complicated manipulation is that of anyons' 
braiding, and its associated non-abelian statistics. While 
in two dimensions the braiding of anyons is defined in 
terms of world lines R(t) that braid one another as time 
evolves, in one dimension - both in the integer m — 1 and 
in the fractional case - a braiding operation requires the 
introduction of tun neling terms between different points 
along the edg d 37 * 38 4 The braiding is then defined in terms 
of trajectories in parameter space, which includes the 
tunneling amplitudes that are introduced to implement 
the braiding. The braiding is topological in the sense that 
it does not depend on the precise details of the trajec- 
tory that implements it , as long as the degeneracy of the 
ground state manifold does not vary throughout the im- 
plementation. Physically, one can imagine realizing such 
operations by changing external gate potentials which de- 
form the shape of the system's edge adiabatically (similar 
to the operations proposed for the Majorana cas d 37 * 39 !). 

In the integer m — 1 case, the interchange of two 
anyons positioned at two neighboring interfaces is carried 
out by subjecting the system to an adiabatically time de- 
pendent Hamiltonian in which interfaces are coupled to 
one another. When two or three interfaces are coupled 
to one another, the degeneracy of the ground state does 
not depend on the precise value of the couplings, as long 
as they do not all vanish at once. Consequently, one may 
"copy" an anyon a onto an anyon c by starting with a 
situation where corresponding interfaces b and c are tun- 
nel coupled, and then turning on a coupling between a 
and b while simultaneously turning off the coupling of b 
to c. Three consecutive "copying" processes then lead to 
an interchange, and the resulting interchanges generate 
a non-abelian representation of the braid group. 

In the integer case only electrons may tunnel between 



two interfaces, thus allowing us to characterize the tun- 
neling term by one tunneling amplitude. In contrast, in 
the fractional case more types of tunneling processes are 
possible, corresponding to the tunneling of any number 
of quasi-particles of charges e/m and spin ±l/m. To 
define the effective Hamiltonian coupling two interfaces, 
we need to specify the amplitudes for all these distinct 
processes. As one may expect, if only electrons are al- 
lowed to tunnel between the interfaces (as may be the 
case if the tunneling is constrained to take place through 
the vacuum), the m = 1 case is reproduced. When sin- 
gle quasi-particles of one spin direction are allowed to 
tunnel (which is the natural case for the FQHE realiza- 
tion of our model), tunnel coupling between either two 
or three interfaces reduces the degeneracy of the ground 
state by a factor of 2m. This case then opens the way 
for interchanges of the positions of anyons by the same 
method envisioned for the integer case. We analyze these 
interchanges in detail below. 

Our analysis of the unitary transformations that corre- 
spond to braiding schemes goes follows different routes. 
In the first, detailed in Section [VI] we explicitly calculate 
these transformation for a particular case of anyons inter- 
change. In the second, detailed in Section [VII| we utilize 
general properties of anyons to find all non-abelian repre- 
sentations of the braid group that satisfy conditions that 
we impose, which are natural to expect from the system 
we analyze. Both routes indeed converge to the same 
result. While the details of the calculations are given in 
the following sections, here we discuss their results. 

To consider braiding, we imagine that two anyons at 
the two ends of the i'th superconducting region are in- 
terchanged. For the m — 1 case the interchange of two 
Majorana fermions correspond to the transformation 
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71 



1 ± i exp(iirQi 



(3) 



This transformation may be written as exp[i^(Qi — k) 2 ] 
with k = 0, 1 corresponding to the ± sign in (|3]), or 
as ^= (1 ± 7172), with 71,72 the two localized Majorana 
modes at the two ends of the superconducting region. Its 
square is the parity of the charge in the superconducting 
region, and its fourth power is unity. Note that in two 
dimensions, the two signs in ^ correspond to anyons 
exchange in clockwise and anti-clockwise sense. In con- 
trast, in one dimension the two signs may be realized 
by different choices of tunneling amplitudes, and are not 
necessarily associated with a geometric notion. Consis- 
tent with the topological nature of the transformation, 
a trajectory that leads to one sign in ^ cannot be de- 
formed into a trajectory that corresponds to a different 
sign, without passing through a trajectory in which the 
degeneracy of the ground state varies during the execu- 
tion of the braiding. 

Guided by this familiar example, we expect that at the 
fractional case the unitary transformation corresponding 
to this interchange will depend only on e 47r< ^ i . We expect 
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to be able to write it as 



2m- 1 

U(Qi)= ^ aj exp (iirjQ; 



(4) 



with some complex coefficients aj, i.e., to be periodic in 
Qi, with the period being 2. We expect the values of aj 
to depend on the type of tunneling amplitudes that are 
used to implement the braiding. 

In our analysis, we find a more compact, yet equivalent, 
form for the transformation U. which is 



U(Qi 



(3) 



The value of a depends on the type of particle which 
tunnels during the implementation of the braiding, while 
the value of k depends on the value of the tunneling am- 

2 

plitudes. For an electron tunneling, a = Just as for 
the m = 1 case, for this value of a the unitary transfor- 
mation |5| has two possible eigenvalues, U 4 — 1, and it 
is periodic in k with a period of 2. For braiding carried 
out by tunneling single quasi-particles we find a = ~ . In 
this case {/ = 1, and U is periodic in k with a period 
of 2m. 

Just as in the m — 1 case, trajectories in parameter 
space that differ by their value of k are separated by tra- 
jectories that involve a variation in the degeneracy of the 
ground state. We note that up to an unimportant abelian 
phase, the unitary transformation Q5J) may be thought of 
as composed of a transformation e la7r ^ !i that results from 
an interchange of anyons, multiplied by a transformation 
e m Q ik that results from a vortex encircling the i'th 
superconducting region 2ak/m times. 

Non-abelian statisti cs is the cornerstone of topological 



quantum computation^ 1 ^ 1 , due to possibility it opens for 
the implementation of unitary transformations that are 
topologically protected from decoherence and noise. It is 
then natural to examine whether the non-abelian anyons 
that we study allow for universal quantum computation, 
that is, whether any unitary transformation within the 
ground state subspace may be approximated by topo- 
logical manipulations of the anyoniP4 We find that, at 
least for unitary time evolution (i.e., processes that do 
not involve measurements) the answer to this question is 
negative, as it is for the integer case. 



III. EDGE MODEL 

The edge states of a FTI a re described by a hydrody- 
namic bosonized theorjBSEE] The edge effective Hamil- 
tonian is written as 



H=™ldx 
2tt 



k(x) (d x <py 



K(x 



(d x ey 



Here, u is the edge mode velocity, <p, 9 are bosonic 
fields satisfying the commutation relation \<f> (x) , 9 (x')] = 
^0 (x' — x) where is the Heaviside step function, and 
gs (x), gF (x) describe position-dependent proximity cou- 
plings to a SC and a FM, which we take to be constant 
in the SC/FM regions and zero elsewhere, respectively. 
The magnetization of the FM is taken to be in the x 
direction. K (x) is a space-dependent Luttinger param- 
eter, originating from interactions between electrons of 
opposite spins. The charge and spin densities are given 
by p = d x 9/ir and s z — d x <fi/ir, respectively (where the 
spin is measured in units of the electron spin H/2). A 
right or left moving electron is described by the opera- 
tors V± = e m ^ ±e \ 

Crucially for the arguments below, we will assume that 
the entire edge is gapped by the proximity to the SC and 
FM, except (possibly) the SC/FM interface. This can be 
achieved, in principle, by making the proximity coupling 
to the SC and FM sufficiently strong. 



IV. GROUND STATE DEGENERACY OF DISK 
WITH 2N SEGMENTS 

We consider a disk with 2N FM/SC interfaces on its 
boundary (illustrated in Fig. [ljj, for N = 2). In order to 
determine the dimension of the ground state manifold, 
we construct a set of commuting operators, which can 
be used to characterize the ground states. Consider the 
operators: e i7rQj = e l ^ i+1 ~ ei \ j = 1, . . . , N, where 9j is 
a 9 field evaluated at an arbitrary point near the middle 
of the jth FM region. The origin (x — 0) is chosen to lie 
within the first FM region (see Fig. [T^). The operator 
6>jv+i is located within this region, to the left of the origin 
(x < 0), while 9\ is to the right of the origin (x > 0). 
The fields 9, (j> satisfy the boundary conditions e l7r< 2 tot = 

e i[6(L-)-6(0+)} and e „5 tot = e i[*(L-)-H0 + )] t where L is 

the perimeter of the system, and Qtot, S'tot are the total 
charge and spin on the edge, respectively. 

Since we are in the gapped phase of the sine-Gordon 
model of Eq. (JsJ) , we expect in the thermodynamic limit 
(where the size of all of the segments becomes large) that 
the 9 field is essentially pinned to the minima of the cosine 
potential in the FM regions. (Similar considerations hold 
for the (j> fields in the SC regions.) In other words, the 
9 — > 9 + 7r/m symmetry is spontaneously broken. In this 
phase, correlations of the fluctuations of 9 decay exponen- 
tially on length scales larger than the correlation length 
£ ~ u/Ap, where is the gap in the FM regions (see 
Appendix [A| for an analysis of the gapped phase). There- 
fore, one can construct approximate ground states which 

= Xj f 0, 



are characterized by (e*^ +1 tijS> ) — (e 17T ® j 
where Xj 



'7.r [g s (x) cos (2m4>) + gF (x) cos (2m9)] . (6) 



A|e™ ?J , where qj € {0, ...,2m— 1} can 
be chosen independently for each FM domain. The en- 
ergy splitting between these ground states is suppressed 
in the thermodynamic limit as e~ R ^, where R is the 
length each region, as discussed below and in Appendix 

El 
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In addition, e MStot commutes both with the Hamilton- 
ain and with e l ^^ j . Therefore the ground states can be 
chosen to be eigenstates of e l7rStot , with eigenvalues e l ™ s , 
s e {0, . . . , 2m — 1}. We label the approximate ground 
states as | {q} ; s) = \q\, . . . , q^; s), where | {q} ; s) satis- 
fies that ({q} ; s\e™ Q i \ {q} ; s) = |A| e&«* . 

For a large but finite system, the |{<?};s) states are 
not exactly degenerate. There are two effects that lift 
the degeneracy between them: intra-segment instanton 
tunneling events between states with different {q}, and 
inter-segment "Josephson" couplings which make the en- 
ergy dependent on the values of {q}. However, both of 
these effects are suppressed exponentially as e~ R ^, as 
they are associated with an action which grows linearly 
with the system size. Therefore, we argue that |{g};s) 
are approximately degenerate, up to exponentially small 
corrections, for any choice of the set {q}, s. 

Similarly, one can define a set of "dual" operators 
e i7rS 3 = e i(^-0 3 _!) 5 j _ 2,...,N, and e l7rSl = 

e twStot Yl i=2 e~ J7rSi . Although the SC regions are in the 
gapped phase, and the fields <f>j are pinned near the min- 
ima of the corresponding cosine potentials, note that the 
approximate ground states | {q} ; s) cannot be further 
distinguished by the expectation values of the operators 
e l7rSj . In fact, these states satisfy ({q}; s|e l7rSj \{q}', s) — > 
in the thermodynamic limit. That is because the oper- 
ators e MrS -» and e i7r< ^ satisfy the commutation relations 

e iTxS iR iTxQj _ e '-^{Si, j + i-Si,j) e iirQj e iirSi ^ ^ 

which can be verified by using the commutation relation 
of the </> and 6 fields. In the state |{g};s), the value of 
giTrQj j g a pp roxuna t e iy localized near e l ™ qj . Applying 
the operator e l7rSj to this state shifts e™Q j to e i7r ^ J+ m), 
as can be seen from Eq. fl7J. This shift implies that 
the overlap of the states \{qj] s) and e™ Sj \ {q}', s) decays 
exponentially with the system size. 

Overall, there are (2m) N+1 distinct approximate eigen- 
state \{q}; s), corresponding to the 2m allowed values of 
charges q^ of each individual SC segment, and the total 
spin s, which can also take 2m values. Not all of these 
states, however, are physical. Labelling the total charge 

by an integer q — X^fLi Qji we see fr° m Eq. (JIJ that s and 
q must be either both even or both odd, corresponding to 
a total even or odd number of fractional quasi-particles 
in the bulk of the system. Due to this constraint, the 
number of physical states is only |(2m) JV+1 . 

In a given sector with a fixed total charge and total 
spin, there are N gs = (2m) JV_1 ground states. For m = 1, 
we get N gs = 2 for each parity sector, as expected 
for 2N Maiorana states located at each of the FM/SC 
interfaces^. 

The ground state degeneracy in the fractional case sug- 
gests that each interface can be thought of as an anyon 
whose quantum dimension is y/2m. This is reminiscent of 
recently proposed models in which dislocations in abelian 
topological phases carry anyons with quantum dimen- 
sions which are square roots of integer d " [ 



V. INTERFACE OPERATORS 

We now turn to define physical operators that act on 
the low-energy subspace. These operators are analogous 
to the Majorana operators in the m — 1 case, in the sense 
that they can be used to express any physical observable 
in the low-energy subspace. They will be useful when 
we discuss topological manipulations of the low-energy 
subspace in the next section. 

We define the unitary operators e lc ^ i and e lv< ^ j such 
that 

e t7r ^^\ qi ,... 7 q N ,s) = e-^-\q x ,...,q N ,s), (8) 



e'^lgi, ...,q N ,s) = \qt, . . . ,q 3 ■, + 1, . . . , q N , s). (9) 

giTrQj j g a diagonal operator in the \{q}, s) basis, whereas 
e l 4>j shifts <7j by one. These operators can be thought 
of as projections of the "microscopic" operators e l ^ j and 
giTrQj ^ introduced in the previous section, onto the low- 
energy subspace. In addition, we define the operator T s 
that shifts the total spin of the system: 

f s \qx,...,q N ,s) = |g 1 ,...,g JV ,s + l}. (10) 

The operators ( |9|10[ ) will not be useful to us, since they 
cannot be constructed by projecting any combination of 
edge quasi-particle operators onto the low energy sub- 
space. To see this, note that they add a charge of 1/m 
and zero spin or spin 1/m with no charge. As a result, 
they violate the constraint between the total spin and 
charge, Eq. (JTJ. However, these operators can be used to 
construct the combinations 

3 

»=i 

X2 J+ i,a = e l ^(T s rl[e^, (11) 

where a = ±1. These combinations, which will be used 
below, correspond to projections of local quasiparticle 
operators onto the low energy manifold. Indeed, the op- 
erators Xj.a {p = il) carry a charge of 1/m and a spin 
of ±l/m (as can be verified by their commutation re- 
lations with the total charge and total spin operators). 
Therefore, their quantum numbers are identical to those 
of a single fractional quasi-particle with spin up or down. 
Moreover, the commutation relations satisfied by Xi,a- 
and for i < j. 

Xi,aXj,l = e l ™Xj,lXi,*, (12) 
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Figure 2. Braiding process, (a) An FTI disk with six 
SC/FM segments. In stages I, II and III of the braiding pro- 
cess, quasi-particle tunneling (represented by blue solid lines) 
is turned on between the SC/FM interfaces, (b) Represen- 
tation of the braiding procedure, involving interfaces 1, 2, 3 
and 4. In the beginning of each stage, the two interfaces con- 
nected by a solid line are coupled; during that stage, the bond 
represented by a dashed line is adiabatically turned on, and 
simultaneously the solid bond is turned off. By the end of 
stage III, the system returns to the original configuration. 



coincide with those of quasi-particle operators e l( t>(x)±ie{x) 
localized at the SC/FM interfaces (for i = j, [xj,ti Xj,i\ — 
if j is odd, and satisfy Xj,tXj,l = e 2in/m Xj,lXj,t if 3 is 
even). Note that in our convention, for j odd, Xj,a cor- 
responds to the interface between the segments labelled 
by e mSj and e n( ^ j , where for j even, between e"® j and 
gttr^+i^ see Fig. [2] a. 

Therefore, the operators Xj,a correspond to quasi- 
particle creation operators at the SC /FM interfaces, pro- 
jected onto the low-energy subspace. This conclusion is 
further supported by calculating directly the matrix ele- 
ments of the microscopic quasi-particle operator between 
the approximate ground states, in the limit of strong co- 
sine potentials (see AppendixjA]) . This calculation reveals 
that the matrix elements of the quasi-particle operators 
within the low-energy subspace are proportional to those 
of Xj,<r , and that the proportionality constant decays ex- 
ponentially with the distance of the quasi-particle oper- 
ator from the interface. We note that the commutation 



relations of Eq. (12 1 ap pear in a one dimensional lattice 
model of "parafermions'SSMl. 



Table I. Summary of the braiding adiabatic trajectory (shown 
also in Fig. [2)3). There are three stages, a = I, II, III, along 
each of which the parameter \ a varies from to 1. The Ham- 
litonian in each stage is written in the middle column, where 
we use the notation Hij = —tijXj,txl -f + h-c. (tij are complex 
parameters). The right column summarizes the symmetry 
operators which commute with the Hamiltonian throughout 
each stage. 



Stage 




Hamiltonian 


Symmetries 


I 


(1 


— Ai) H12+X1H23 


g!T<23 e i7rS 3 


II 


(1- 


- An) H23+X11H24 




III 


(1- 


Am) H24+Ani-ffi2 





inal form, but the state of the system does not return to 
the initial state. The adiabatic evolution corresponds to 
a unitary matrix acting on the ground state manifold. 

Below, we analyze a braid operation between nearest- 
neighbor interfaces, which we label 3 and 4 (for later con- 
venience). The operation consists of three stages, which 
are described pictorially in Fig. [2]d. It begins by nucle- 
ating a new, small, segment which is flanked by the in- 
terfaces 1 and 2. At the beginning of the first stage, the 
small size of the new segment means that interfaces 1 and 
2 are coupled to each other, and all the other interfaces 
are decoupled. During the first stage, we simultaneously 
bring interface 3 close to 2, while moving 1 away from 
both 2 and 3, such that at the end of the process only 
2 and 3 are coupled to each other, while 1 is decoupled 
from them. In the second stage, interface 4 approaches 
3, and 2 is taken away from 3 and 4. In the final stage, 
we couple 1 to 2 and decouple 4 from 1 and 2, such that 
the Hamiltonian returns to its initial form. In the follow- 
ing, we analyze an explicit Hamiltonian path yielding this 
braid operation, which is summarized in Table ^ Later, 
we shall discuss the conditions under which the result is 
independent of the specific from of the Hamiltonian path 
representing the same Braid operation. 



B. Ground state degeneracy 



VI. TOPOLOGICAL MANIPULATIONS 
A. setup 

The braiding process is facilitated by deforming the 
droplet adiabatically, such that different SC/FM inter- 
faces are brought close to each other at every stage. Prox- 
imity between interfaces essentially couples them, by al- 
lowing quasi-particles to tunnel between them. We shall 
assume that only one spin species can tunnel between 
interfaces. The reason for this assumption will become 
clear in next sections, and we shall explain how it is man- 
ifested in realizations of the model under consideration. 
At the end of the process, the droplet returns to its orig- 



To analyze the braiding process, we first need to show 
that it does not change the ground state degeneracy. We 
consider a disk with a total of N = 3 segments of each 
type. The ground state manifold, without any coupling, 
is (2m) 2 fold degenerate. We define operators #12, #23 
and H24, the Hamiltonians at the beginning of the three 
stages I, II, III. These are given by 

H jk = -tjkXj, t >d,t + h - c - ( 13 ) 

where the tjk are complex amplitudes. 

Consider first the initial Hamiltonian (see Table |l|, 
given by 
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#i2 = -foxa.txt.t + h - c - = ~ 2 l il2 ' cos ( 7r( 3 1 + 

(14) 

Here, tp±2 — arg (£12). It is convenient to work in the ba- 
sis of eigenstates of the operators e I7r( ^ 1 , e l7r ® 2 , e l7r< ^ 3 , and 
e" Stot , which we label by I91, 92, 93! s)- The total charge 
and spin are conserved, and we may set . q^- = and 

s = 0. Then, a state in the (2m) 2 -dimensional low-energy 
subspace can be labelled as | (72,93), where 91 is fixed to 
qi = — 92 — 93 ■ The initial Hamiltonian (14 1 is diagonal in 



this basis, and therefore its eigen-energies can be read off 
easily: E 12 (q 2 ,q 3 ) = -2 |t 12 | cos [-^ (q 2 + 93) + (p 12 \ ■ 
For generic <pi2, there are 2m ground states. Since this 
FM segment is nucleated inside a SC region, its total spin 
is zero, and the ground states are 



l**(?3)) = 1 92 = -93,93), 



(15) 



labelled by a single index 93 = 0, ... , 2m — 1. The resid- 
ual 2m-fold ground state degeneracy can be understood 
as a result of the symmetries of the Hamiltonian. From 
Eq. (141 e i7T ® 3 and e"^ 3 commute with H 12 . The com- 



mutation relations between e I7r< ^ 3 and e I7rSs ensures that 
the ground state is (at least) 2m-fold degenerate by the 
eignevalues of e l7r ® 3 . 

Similar considerations can be applied in order to find 
the ground state degeneracy throughout the braiding op- 
eration. The operator e Z7T< ^ 3 always commutes with the 
Hamiltonian, at any stage. This can be seen easily from 
the fact that the segment labelled by e l7r< 3 3 never couples 
to any other segment at any stage (see Fig. [2^). Using 
the definition of the \i<y operators, Eq. (Ill, one finds 
that 



#23 = -2|t 2 3|cos 



I cos (ttS 2 + V23) 



and 



#24 = -2 |*24 I COS 7T ( S 2 + Q 2 I + 1^24 



(16) 



(17) 



In each stage, a = I, II, III, there is a symmetry op- 
erator S Q that commutes with the Hamiltonian, and 
satisfies E a e in ® 3 = e -i £e^ ! E a . We specify S Q for 
each stage in the right column of Table |TJ and the afore- 
mentioned relation can be verified using Eq. This 
combination of symmetries dictates that every state is at 
least 2m fold degenerate, where each degenerate subspace 
can be labelled by q 3 . Assuming that the special values 
^12,9523 = 7r(2i + l)/(2m) and (/J24 = wl/m (I integer) are 
avoided, the ground state is exactly 2m-fold degenerate 
throughout the braiding process (the special values for 
the (fiij give an additional two fold degeneracy). Note 
that these conclusions hold for any trajectory in Hamil- 
tonian space, as long as the appropriate symmetries are 
maintained in each stage of the evolution, and the acci- 
dental degeneracies are avoided. 



C. Braid matrices from Berry's phases 

The evolution operator corresponding to the braid op- 
eration can thus be represented as a block-diagonal uni- 
tary matrix, in which each (2m) x (2m) block acts on 
a separate energy subspace. We are now faced with 
the problem of calculating the evolution operator in the 
ground state subspace. Let us denote this operator by 
i/34, corresponding to a braiding operation of interfaces 
3 and 4. The calculation of ^34 can be done analytically 
by using the symmetry properties of Hamiltonian at each 
stage of the evolution. 

We begin by observing that, since e Z7r( ^ 3 always com- 
mutes with the Hamiltonian, U34 and the evolution op- 
erators for each stage are diagonal in the the basis of 
e *7T<23 eigenstates. In every stage, the adiabatic evolu- 
tion maps e Z7rC " eigenstates between the initial and final 
ground state manifolds while preserving the eigenvalue 
q 3 , and multiplies by a phase factor that may depend on 
93. This is explicitly summarized as 



#a|*?(9 3 )> = exp (»7« (93)) |*?(<fe)>. (18) 

Here, U a is the evolution operator of stage a = I, II, III, 
and (93)) are t ne ground states of the initial (fi- 

nal) Hamiltonian in stage a, respectively, which are la- 
belled by their e iTr< ^ 3 eigenvalues. Likewise, j a (q 3 ) are 
the phases accumulated in each of the stages. 

In order to determine 7^ (93), we use the additional 
symmetry operator S Q for each stage, as indicated in 
Table HJ This symmetry commutes with the Hamiltonian, 
and therefore also with the evolution operator for this 
stage [E a ,[/ Q ] = 0. Acting with S a on both sides of 
( 18 1, we get that 



U a X a \H>?(q 3 )) = e^(«)S a |^(9 3 )). (19) 

Furthermore, the relation e I7r< ^ 3 S a = e l ™ S Q e l7r< ^ 3 im- 
plies that the operator E a advances e I7r( ^ 3 by one incre- 
ment, and therefore for both the initial and final stage at 
each stage we have, 



Oft)) = exp (iSfa (9 3 ))|*? (/) (93 + 1)), (20) 

where $frf\ (93) are phases which depends on gauge 
choices for the different eigenstates, to be determined be- 
low. Inserting ( 20 1 into ( 19 ) , we get the recursion relation 



7a (93 + 1) = 7a (93) + 8? (93) " 5? (93) • 



(21) 



Note that while the phase accumulation at each point 
along the path depends on gauge choices, the total Berry 



phase accumulated along a cycle does not. It is conve- 
nient to choose a continuous gauge, for which the total 
Berry's phases are given by 

Us4n(Qs)) = exp(iJ2la(q 3 ))\mQs))- (22) 

a 

A continuous gauge requires {^(ns)) — (ns)} . 

Therefore, the values of the phases 8f,^ depend only on 
three gauge choices. These are the gauge choices eigen- 
states of the Hamiltonians H\2, H23, and H24, which con- 
stitute the initial Hamiltonian at the beginning of stages 
I-III, as well as the final Hamiltonian for stage III. 

Making the necessary gauge choice, allows us to solve 
Eq. (21 1 for 70(^3), yielding the total Berry phase (the 
details of the calculation are given in Appendix [B| 



2m 



(<?3 - kf 



(23) 



The integer k depends on the choice for the phases (fij. 
Recall that the Hamiltonians Hij , Eqs. ( 14 1 , ( 16 - 17 1 , have 
an additional degeneracy for a discrete choice of the (fij . 
Any two choices for the ifiij that can be deformed to each 
other without crossing a degeneracy point yield the same 
k. 

The evolution operator for the braiding path can 
be written explicitly by its application on the eigen- 
states of the Hamiltonian in the beginning of the c ycle , 



U. 



i^te)) = e 



*i((fe)>. Since by Eq. (151, 



the ground states of the initial Hamiltonian satisfy 52 = 
— qs, this can be written in a basis-independent form in 
terms of the operator e l7r( ^ 2 . Loosely speaking, i/34 can 
be written as 



U. 



34 



exp — -— Q2 



(24) 



Alternatively, using the identitjE^ e*a 



3777''" 



1 sr^2rn- 
2m 2^p=0 



2m-l J% [pg-Trj+il 



one can write 



U. 



Ik) 



31 



2m- 1 

E e " 

p=0 



(25) 



In the case m = 1, t/34 reduces to the braiding rule of 
Ising anyons^HUl 

Following a similar procedure, one can construct the 
operator representing the exchange of any pair of neigh- 



boring interfaces: U^j-i 2j 



■-(Sj-k/m)- Mk) 



2j,2j+l 



e 2 (Oj+i-fc/m) _ j n orc | er f or these operations to form 
a representation of the braid group, it is necessary and 
sufficient that they satisfy 



U, 



(fc») M k i) ' 



(N-i|>i), 



(26) 



r(fei) 



J j+l,j+2 u j,j+l 



iS J+9 - ( 27 ) 





Q 



Figure 3. Diagrammatic representation of the Yang- 
Baxter equations (Eq. D2|. Three interfaces 1,2,3 are 
braided in two distinct sequences. The Yang-Baxter equa- 
tions state that the results of these two sequences of braiding 
operations are the same. 



Equation ( 26 1 clearly holds because the spin or charge op- 
erators of non-nearest neighbor segments commute. Us- 



ing ( 25 ) , it is not difficult to show that ( 27 1 holds as well 
(see Appendix [D] Eq. ( 27 1 is depicted in Fig. [3]) . There- 
fore, Ui^+i form a representation of the braid group. In 
that respect, our system exhibits a form of non-abelian 
statistics. By combining a sequence of nearest-neighbor 
exchanges, an exchange operation of arbitrarily far seg- 
ments can be defined. 

In any physical realization, we do not expect to con- 
trol the precise form of the Hamiltonian in each stage. 
It is therefore important to discuss the extent to which 
the result of the braiding process depends on the details 
of the Hamiltonian along the path. We argue that the 
braiding is "topological", in the sense that it is, to a large 
degree, independent of these precise details. 

To see this, one needs to note that the braiding unitary 
matrix was derived above without referring to the precise 
adiabatic path in Hamiltonian space. All we used were 
the symmetry properties of the Hamiltonian in each stage 
( Table [I]). These symmetries do not depend on the precise 
details of the intermediate Hamiltonian, but only on the 
overall configuration, e.g., which interfaces are allowed to 
couple in each stage. 

In Appendix |C 1| we state more formally the conditions 
under which the result of the braiding is independent 
of details. Special care must be taken in stage III of 
the braiding, in which quasi-particles of only one spin 
species, e.g. spin up, must be allowed to tunnel between 
interfaces 2 and 4. We elaborate on the significance of 
this requirement and the ways to meet it in the various 
physical realizations in Appendix |C 2| 



VII. 



BRAIDING AND TOPOLOGICAL SPIN OF 
BOUNDARY ANYONS 



In the previous section, we derived the unitary matrix 
representing braid operations by an explicit calculation. 
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In the following, we shall try to shed light on the physical 
picture behind these representations. To do so, we show 
that the results of the previous section can be derived 
almost painlessly, just by assuming that the representa- 
tion of the braids have properties which are analogous to 
those of anyons in two dimensions. The first and most ba- 
sic assumption, is very natural: there exists a topological 
operation in the system which corresponds to a braid of 
two interfaces, in that the unitary matrices representing 
this operation obey the Yang-Baxter equation. 

The operations we consider braid two neighboring in- 
terfaces, but do not change the total charge (spin) in the 
segment between them. This results from the general 
form of the braid operations - to exchange two interfaces 
flanking a SC (FM) segment, we use couplings to an aux- 
iliary segment, of the same type. Therefore, charge (spin) 
can only be exchanged with the auxiliary segment. Since 
the auxiliary segment has zero charge (spin) at the be- 
ginning and end of the operation, the charge of the main 
segment cannot change by the operation. Indeed, this 
can be seen explicitly in the analysis presented in the 
previous section. As a result, the unitary matrix repre- 
senting the braid operation is diagonal with respect to 
the charge (spin) of the segment. 

The derivation now proceeds by considering a property 
of anyons called the topological spin (TS). In two dimen- 
sions, the topological spin gives the phase acquired by a 
2ir rotation of an anyon. For fermions and bosons, the 
topological spin is the familiar —1 and +1 respectively 
(corresponding to half-odd or integer spins). There is a 
close connection between the braid matrix for anyons and 
their topological spin. In two dimensions, thes e relations 
have been considered by various authors 9 4? . The sys- 
tem under consideration is one dimensional, and there- 
fore seemingly does not allow a 2tt "rotation" of a particle. 
However, as we shall explain below, the TS of a particle 
can be defined in our system using the relations of the TS 
to the braid matrix. We shall then see how to use these 
relations to derive the possible unitary representations of 
the braid operations in the system at point. 

In our one dimensional system, we consider the TS 
of two different kinds of objects (particles) - interfaces, 
which we denote by X, and the charge (or spin) of a 
segment, which we shall label by q = 0,1, ..2m — 1. In 
what follows, we need to know how to compose, or fuse 
different objects in our system. As we saw above, two 
interfaces yield a quantum number exp(i7rQ) (exp(z7rS')) 
which is the total charge (spin) in the segment between 
them, respectively. Suppose we consider two neighbor- 
ing SC segments with quantum numbers exp(inQi) and 
exp(i7rQ2), and we "fuse" them by shrinking the FM re- 
gion which lies between them. This results in tunneling 
of fractional quasi-particles between the two SC regions, 
and energetically favors a specific value for exp(i7r5) in 
the FM region. The two SC segments are for all pur- 
poses one, where clearly, in the absence of other cou- 
plings, exp (i7r(Qi + Qi)) remains a good quantum num- 



ber. This therefore suggests the following fusion rules 



YxY = 0+l + ... + 2m-l 
Qi x Q2 = (<7i + 92) mod 2m 



(28) 



We note that the labelling qi = 0, 1..., 2m— 1 does not de- 
pend on the gauge choices in the definition of the opera- 



tors exp(i ttQ). Equation (28 1 suggests that the labelling 
can be defined by the addition law for charges, in which 
each type of charge plays a different role. Indeed, this 
addition rule has a measurable physical content which 
does not depend on any gauge choices. 

In two dimensional theories of anyons, it is convenient 
to think about particles moving in the two dimensional 
plane, and consider topological properties of their world 
lines (such as braiding). In this paper, we have de- 
fined braiding by considering trajectories in Hamiltonian 
space. In the following, we represent these Hamiltonian 
trajectories as world-lines of the respective "particles" in- 
volved, keeping in mind that they do not correspond to 
motion of objects in real space. 

We are now ready to define the TS in our system. In 
short, the TS of a particle is a phase factor associated to 
the world line appearing in FigQa). For interfaces, it is 
concretely defined by the phase acquired by the system 
by the following sequence of operations, as illustrated in 
Fig |4^e): (i) nucleation of a segment to the right (by 
convention) of the interface X\ (note that the notation 
Xi corresponds to particle X at coordinate r%). The 
total spin or charge of this segment is zero (the nucle- 
ation does not add total charge to the system). The 
couplings between X2 and A3 flanking the new segment 
is taken to zero, increasing the ground state degeneracy 
by a factor of 2m. (ii) A right handed braid operation 
is performed between Xi and A 2 . (iii) The total charge 
q of the segment between X2 and A3 is measured, and 
we consider (post-select) only the outcomes correspond- 
ing to zero charge. Therefore, the system ends up in 
the same state (no charges have been changed anywhere 
in the system), up to a phase factor. Importantly, this 
phase factor does not depend on the state of the system, 
since the operation does not change the total charge in 
the segment of X\ and A2 (see Appendix [E] for a more 
detailed discussion). We can therefore define this phase 
factor as Ox, the topological spin of particle type X. 

In order to define q , the topological spin of a charge 
q, we first need to define the operation corresponding to 
an exchange of two charges. Consider the sequence of 4 
right handed exchanges of the interfaces, as in Fig. |4](c). 
The figure suggests that this sequence should yield an 
exchange of the fusion charges q\ and qi of the two pairs 
of X particles, as would indeed be the case for anyons 
in two dimensions (see Appendix. [E] for more details). 
Our second assumption is that this is indeed the case. 
Since by Eq. ( 28 1 there is only one fusion channel for 



the gi's, the state is multiplied by a phase factor which 
depends only on q\ and qi - the charges qi are abelian. 
It is straightforward to check that exchanges of charges 
satisfy the Yang-Baxter equation. 
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The operations defining 9 q are illustrated in Figs |4^b) 
and (f). We consider, say, a SC segment in an eigenstate 
of exp(i7rQ 1 ) to which we associate a particle type q 1 . In 
step (i), we nucleate two segments to the right of q\. This 
step is actually done in two substeps: first, we nucleate 
a SC segment to the right of the segment qi, and then 
we nucleate a FM segment which separates this segment 
into two segments, with charges exp (iir (Q 2 + Q3)) = 1. 
Note that the spin of the middle FM segment is also 
zero, exp(i7rS f 2) = 1. We now perform a braid between 
the segments labelled qi and q 2 . Next, we measure the 
charges exp (in (Q 2 + Q3)) and exp(wr,f?), post-selecting 
their values to be equal to 1. The charge q± in the orig- 
inal segment is therefore unchanged by this sequence of 
operations, as well as any other charges used to label 
the original state of the system. As before, the state of 
the system acquires a phase, which only depends on the 
charge q\. 

We now use an important relation between the TS of a 
composite to a double braid of its two components. This 
is a result of an equality between Fig|4](b) and (d), which 
can be derived by using the Yang-Baxter equation and 
the definition for the TS. The equality between Fig.[4](b) 
and (d), with both incoming and outgoing lines labelled 
by X. yields 



■U\q) 



(29) 



where U(q) corresponds to the unitary matrix represent- 
ing an exchange of two interfaces whose fusion charge is 
q. We shall not attempt to calculate Ox- However, by 
calculating the topological spin 9 q of the composite, we 
can get, using Eq. (29), the square of the sought after 



braid matrix, up to a global phase. 

In two dimensional theories of anyons, the TS of a com- 
posite of two particles is equal to the TS of the particle 
they fuse into. This can be understood by noting that 
the equality between Fig[4](b) and (d) guarantees that in 
two dimensions, the TS of the composite includes both 
the intrinsic spin of the two particles and their relative 
angular momentum (see Appendix [E] for more details) . 
We now assume the same holds in our one dimensional 
system - the TS of the qi, q 2 composite is equal to the 
TS of combined q = q\ + q 2 mod 2m charge. The TS 
of a composite of two charges q\ and q 2 is just the pro- 
cess described in Fig. Qb), with the outgoing lines are 
labelled with q\ and q 2 . 

Next, we note that for consistency, the charge q = 
must correspond to trivial TS, #0 = 1- We take the TS 
of the elementary charge q — 1 as a parameter, 9i = e la , 
to be determined later. Importantly, note that Fig. Qa) 
implies that an exchange of two segments with q = 1 
leads to the phase factor B\ . Using the composition rule 
in Eq. ( 28 ) , 6 q should be equal to the TS of a composite 



of q unit charges, a process where we encounter q ex- 
changes of these elementary charges. This gives for the 
topological spin 




Figure 4. Topological Spin. The process illustrated in (a) 
defines the TS, where the cup U corresponds to creation of 
a particle -anti particle pair (two X's fusing to q = or two 
charges with 51+92 = 0), and the cap f] corresponds to 
projection on zero total charge, (b) TS of a composite ob- 
ject. When both lines are labelled by X and fuse to charge 
q , the phase acquired is 9 q . The 4-fold crossing, which is 
magnified in (c), is assumed to result in an exchange of two 
q's. When the lines are labelled by qi and 92, the phase ac- 
quired is 6 qi+qo mo d 2m- Importantly, the TS of a composite 
is equal to the phase accumulated in the process appearing in 
(d). This equality results from using the Yang-Baxter equa- 
tion and the definition of the TS. (e) The evolution of the 
different segments in the process defining 8x- The green ar- 
row corresponds to braiding, while the "M" corresponds to 
projection on zero charge, (f) The process defining 9 q . The 
braids are ordered by color, green, blue and finally purple. 
The projections are on zero spin for the magnet segment in- 
tersecting the two SC segments, and zero total charge for 
these two segments. 



The structure above appears since fusing 2m ele- 
mentary particles results in the trivial q = charge, and 
requires that 



up 
2tt? 



pez. 



Taking the square root of Eq. ( 30 ) we arrive at 



U(q) = (-1) 



/(<?) 



exp 



(31) 



(32) 



where f(q) is an arbitrary, integer valued function. To 
determine a and f(q) we appeal to the Yang-Baxter equa- 
tion, Eq. (27 1. We find numerically that solutions are 



possible only for even p's in Eq. (311, where sign func- 
tion (— l)-^ 9 ) can be of the form e 47r ("i9 +"2) with integer 
m,ri2- Therefore, the extra (—1) signs can always be ab- 
sorbed in the definition of a, up to an overall sign. The 
final result is therefore 



0q mod 2m = CXp (iaq 2 



(30) 



U(q) =exp(i^VY 
V 2m / 



n € Z. 



(33) 
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As noted in Sec. [TTJ the value n = 1 is realized by 
quasi-particle tunnelling of a single spin species, while 
n = m is realized by electron tunnelling. The other 
representations of the braid group, which are given by 
U^ k \q) = exp (z^-(<7 — fc) 2 ) can be obtained from the 
above considerations by adding more particles types. In- 
deed, we see that 



(34) 



is, up to a global phase, a topological operation which 
is equivalent to taking a vortex n times around a SC 
segment, or changing the orientation of the FM direction 
n times in the x — y plane by 2ir. 



VIII. QUANTUM INFORMATION 
PROCESSING 



In order to assess the suitability of our system for topo- 
logically protected quantum computation, it is instruc- 
tive to examine the structure of the resulting non-abelian 
theory. Below, we show that the representation of the 
braid group realized in our system is a direct product of 
an Ising anyonic theory times a novel representation of di- 
mension m. We argue that braiding operations alone are 
not sufficient to realize universal topologically protected 
quantum computation, in agreement with the general ar- 
gument for models with anyons of quantum dimension 
which is a square root of an integer^. 

The unitary matrix that describes the braiding oper- 
ation of two interfaces at the ends of a superconduct- 
ing segment, Eq. (24 1, depends on the charge in the 
superconducting segment. The charge can be written 
as Q = q/m where q can be uniquely expressed as 
q = mq a + 2q v , with q a — 0, 1 and q v — 0, . . . , m — 1. 



Inserting this expression into Eq. (24 1, and assuming for 
simplicity k = 0, we get 



U 



e 2 



(35) 



Therefore, we see that if we write the Hilbert space % 
as a tensor product % a ® H v , such that the states are 
written as \q) — \q a ) ® \q v ), the braiding matrix decom- 
poses into a tensor product U a <g> U v . Here, U a and U v 
are 2x2 and m x m matrices, given by the first and 
second terms on the right hand side of Eq. ( [35| , respec- 
tively. A similar decomposition holds for a braid op- 
eration acting on a ferromagnetic segment. In this re- 
spect, we see that the 2m-dimensional representation of 
the braid group given by Eq. (24 1 is reducible: it de- 



composes into a two-dimensional representation, which 
is nothing but the representation formed by Ising anyons 
cr, times an m-dimensional representation corresponding 
to braiding of a non-Ising object v. 

This decomposition gives insights into the class of uni- 
tary transformations that can be realized using braid- 
ing of interfaces, and hence their suitability for quantum 



computation. We now argue that by using braiding op- 
erations alone, the system studied in this paper does not 
allow for topologically protected universal quantum com- 
putation. Ising anyons are known not to provide uni- 
versality for quantum computatio n 129 ! 30 ^. Due to the 
tensor product structure of the topologically protected 
operations, it is sufficient to consider the fractional part, 
corresponding to Ti. v . The braid representations acting 
within this subspace preserves^ a generalization of the 
Pauli group to q-dits of dimension m. Therefore, the 
braid operations can be simulated on a classical com- 
puter, and are not universal. 

Conceptually, universality could be achieved by adding 
an entangling operation between the Ising and the frac- 
tional parts. A braiding operation would then produce 
an effective phase gate that would provide the missing 
ingredient to make the Ising part universal. However, at 
present we do not know whether it is possible to realize 
such an operation in a topologically protected way. More- 
over, topologically protected measurements of charges^! 
(and spins) cannot achieve such entanglement, since sim- 
ilarly to the braiding, they can also be shown to be of a 
tensor product form. 



IX. CONCLUDING REMARKS 

In this work, we have described a physical route for 
utilizing proximity coupling to superconductors in order 
to realize a species of non-abelian anyons, which goes 
beyond the Majorana fermion paradigm. The essential 
ingredients of the proposed system are a pair of counter- 
propagating edge modes of a Laughlin fractional quan- 
tum Hall state, proximity-coupled to a superconductor. 
As we saw, there are several possible realizations of such 
a system. One could start from a "fractional topolog- 
ical insulator" whose edges are coupled to an array of 
superconductors and ferromagnets. In the absence of 
any known realization of a fractional topological insu- 
lator phase (as of today), one could get by starting from 
"ordinary" Laughlin fractional quantum Hall state whose 
edges are coupled to a superconductor. The fractional 
quantum Hall state in graphene might be a promising 
candidate for realizing such systems, since the magnetic 
fields needed for observing it are much lower than the 
fields needed in semiconductor heterostructure devices. 

An experimentally accessible signature of the fraction- 
alized Majorana modes is a fractional Josephson effect, 
which should exhibit a component of 4to7t periodicity 
(analogously to the 4tt periodicity predicted for topo- 
logical superconductors with Majorana edge modesM^ 
In addition, it might be possible to observe topological 
pumping of fractional charge by controlling the relative 
phase of the superconducting regions. 

More broadly speaking, the system we describe here is 
an example of how gapping out the edge state of a frac- 
tionalized two-dimensional phase can realize a topological 
phase which supports new types of non-abelian particles, 



13 



not present in the original two-dimensional theory. In 
our example, the underlying Laughlin fractional quan- 
tum Hall state supports quasiparticles with a fractional 
charge and fractionalized abelian statistics; the resulting 
gapped theory on the edge, however, realizes non-abelian 
quasiparticles. Moreover, the resulting non-abelian the- 
ory on the edge is shown to go beyond the well-known 
Majorana (Ising) framework. 

This may seem contradictory to the general 
axgumenta^HSS] indicating that gapped one-dimensional 
systems with no symmetry other than fermion parity 
conservation support only two distinct topological 
phases, a trivial phase and a non-trivial phase, with 
an odd number of Majorana modes at the interface 
between them. The reason our system avoids this ex- 
haustive classification is that it is not, strictly speaking, 
one-dimensional; the edge states of fractional quantum 
Hall states can never be realized as degrees of freedom 
of an isolated one-dimensional system. This is reflected, 
for example, in the fact that the theory contains 
"local" (from the edge perspective) operators which 
satisfy fractional statistics, which is not possible in any 
one-dimensional system made of fermions and bosons. 

It would be interesting to pursue this idea further, 
by examining gapped states which are realized by gap- 
ping out edge modes of topological phases. This may 
serve as a route to discovering new classes of topological 
phases with non-abelian excitations. For example, more 
complicated^ quantum Hall states or higher dimensional 
fractional topological insulators^HUl mav b e interesting 
candidates for such investigations. 



Appendix A: Appendix: Matrix elements of the 
quasi-particle operator 

In this Appendix, we describe an explicit calculation of 
the matrix element of a quasi-particle operator between 
different states in the ground state manifold. We show 
that the matrix clement is finite if the quasi-particle is 
located sufficiently close to an interface between a su- 
perconducting and a ferromagnetic segment. The matrix 
element decays exponentially with the distance from the 
interface. 



1. Model 

Let us consider a system composed of one supercon- 
ducting segment, extending from x ~ to x — L, be- 
tween two long ferromagnetic segments at x < and 
x > L. For simplicity, we assume that the gap in the 
ferromagnetic segments is very large, such that charge 
fluctuations are completely quenched outside the super- 
conductor. The Hamiltonian for < x < L is 



// = J d.r 



{^u[{d x <j)f + {d x ey 



gs cos (2m<f>)\ 
(Al) 



supplemented by the boundary condition 



d x <t>{x = Q,L) = 0, 



(A2) 
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which accounts for the fact that the current at the edges 
of the superconductor is identically zero, due to the large 
gap in the ferromagnetic regions. We are assuming that 
the coupling gs is large enough such that the field (f) is 
pinned to the vicinity of the minima of the cosine poten- 
tial, <f> ~ ^l, where I is an integer. Deep in the super- 
conducting phase, one can expand the cosine potential 
up to second order around one of the minima, obtaining 
the effective Hamiltonian 









m 


— u 




2tt 



where g = (2m) gs- 



[d^f + {d x of 









m ) 



The Hamiltonian (A3) is quadratic, and can be diago- 



nalized using the following mode expansion: 



.(x) = — 1 + 

771 



V k=l 



cos 



7rfc 



(a k - 4) > 
(A4) 
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'(*) = 0(O)- 



mL 



k=l 



sm 



irk 



(o* + 4) 



manifold, which wc index by the average values of and 
(j) on either side of the x = interface. A diagonal matrix 
element has the form 



(A5) 



Here, we have introduced the ladder operators a k , k 



1,2, ... , satisfying 



5k,k' and [a fe ,a fc /] = 0. </> 



and h are the average phase and the charge of the super- 
conducting segment. These variables are canonical conju- 



gates, satisfying 



& ,n 



Note that [9 (0) , H, 



off 



0, 



therefore 9 (0) can be replaced by a c-number: 
0(0) = ~P, 



(A6) 



with integer p, where we have assumed that these values 
minimize the (infinite) cosine potential on the ferromag- 
netic side x < 0. Using Eq. ( A4 1 , ( A5 ) , one can reproduce 
the commutation relation [</> (x) , 6 (x 1 )] = i — Q (x' — x). 
Inserting the mode expansions into the Hamiltonian 



( A3 ) , we get 



a\a k + a k a\ 



ttu o qL 79 / Tiuk qL \ 

k=l 

This Hamiltonian is diagonalized by a Bogoliubov trans- 
formation of the form 

a k = a k b k + I3 k bl, (A8) 

where a k = ^(f + l), ft = -l), and 

£ fc = ^/A* - B*, expressed via A k = ^ + ^ and 
= ^-r. The n, <^o part of i? c ff is diagonalized by 
introducing ladder operators 77, such that 



\ TTU 

( graL 
n = 

\ TTU 



-1/4 



V + V 



it 



2\ 1/4 « ~+ 
- \ ^ _ T ^T 



The diagonal form of H e g (up to constants) is 



sc 



k=l 



(A9) 



(A10) 



where _B = yfirgu/r 



i[4>{x)+6{x)\ 



VP. 



l,p) 



(All) 



where \i^i, p ) — \itl/m, itp/m) is a ground state in which 
4> (x > 0) and 6> (x < 0) are localized near — I and ^p, re- 
spectively. Note that these two variables commute, and 
therefore they can be localized simultaneously. To eval- 
uate A^p (x), we use the identity 



3 dv = e <o) + i«d 2 )-<6) 2 ) 



(A12) 



valid for any operator O which is at most linear in 
creation and annihilation operators. Substituting O — 
i[<f)(x) + 9 (x)], the expectation values in the exponent 
can be computed using the mode expansions (A4|, (A5l. 



The computaion is lengthy but straightforward, giving 



(x) 



(A13) 



where 



(mg\ -V 2 



2L \ ttu / 



k=i 



TTU COS' 



[fx] 



■Kuk ^ 
L y 



(A14) 



TTvgu 



F e {x) = [ — J — — — 

V TTU I IL 



£ 

k=l 



TTUk 

L 



ink 



sm 



' irk - 



/ TTUk 

V L 



(A15) 



Here, we have introduced exponential damping factors 
of the form , where a is a short-distance cutoff, to 
suppress ultraviolet singularities. (x) — (cf> 2 (x)) — 
(4> (x)} 2 , and similarly for Fg. (Expectation values of the 
form (<f> (x) 9 (x)) vanish.) 

We now analyze the asymptotic behavior of F^ and Fg 
for £ <C x <C L. where we have defined the correlation 
length as £ = ^Jmuj (rrg). In the limit L — >• 00, the 
sums over k in Eq. ( |A14 ),( A15 1 can be replaced by inte- 
grals over q = k/L. Then, the long-distance asymptotic 
behavior of F^ and Fg is easily extracted: 



(A16) 



2. Computation of the matrix elements 

Next, we calculate matrix elements of a quasi-particle 
creation operator between states in the ground state 



Fe{x) 



1 

2mf 



1 a* 

--alog- 



4x 2 



ix log 



a — 2zx x 

a + 2ix , 

(A17) 
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Inserting these expressions into (A13) gives 



Ai, P (x) 



f('+p)- 



(A18) 



Therefore, the diagonal matrix element of the quasi- 
particle operator in the ground state |i/';,p) decays expo- 
nentially with the distance from the interface. In a very 
similar way. one can show that the matrix element of the 
quasi-particle operator between two ground states with 
different {l,p} vanishes in the limit L —¥ oo. It is there- 
fore natural to identify the operators \ i (T introduced in 
Sec. [V] as the projection onto the ground state subspace 
of the quasi-particle operators e*^ x ' ±e ™ acting at the 
interface. 



Appendix B: Calculation of the braid matrix 

In order to complete the calculation of the unitary ma- 
trix corresponding to the braiding operation of two in- 
terfaces, we first need to obtain the ground states of the 
Hamiltonians -H12, H23 and H24, making the necessary 
gauge choices. As in Sec. |VI| we use the basis of eigen- 
states of the operators exp(inQj), j — 1,2,3. We work in 
the sector . exp(iirQj) = 1. A state in this sector can 
be labelled as |<?2,93), where exp(^qj) is the eigenvalue 
of exp(mQj) and q x = — q-i — 33. 

The ground state of H 12 = — 2|t 12 | cos(7rQi + V12) is 
given by 

|$(«s)) = l4 n (<Z 3 )> = \<b = -Q3 + ki,qa), (Bl) 
where the integer k\ is determined by ip\2 according to 



— (fci- -) < tp ia < — 
ml ml 



(B2) 



The ground state is 2m-fold degenerate, corresponding to 
the 2m possible values of 53. Equation (Bl) includes an 



explicit gauge choice for the ground states. Note that for 
= —(ki + |)) the ground state degeneracy increases 
to 4m. We therefore assume that these values of ipi2 are 
avoided. 

The Hamiltonian in the beginning of the second stage, 
H23, can be written in the basis of g 2 eigenstates as 



Importantly, note that when = ir, the ground state 
of H23 is doubly degenerate. These are degeneracy points 
that we assume are avoided in the braiding process. 

The ground state for a particle on a ring with flux <&^g 
is simply a plane wave, 



2m-l 

E e ™ n l92 

n=0 



-q3+n,q 3 ), 



(B5) 

where ku is the closest integer to — (771(^23) /tt- Note that 
again, a gauge choice for the overall phase of the states 
has been made in Eq. (|B5| 



The phases 8\ /fe), defined in Eq. (20 1 of the main 
text, are determined by operating with the symmetry 
operator Si = exp^iS^) on the ground states of the ini- 
tial (H12) or the final (H23) Hamiltonian, Eqs. (Bl I , ( B5 1, 



respectively. In the gauge we have chosen, this gives 

Sj (q 3 ) = 6f (q 3 ) = 0. (B6) 
Therefore, the recursion relation for 70(93), 

7a (q 3 + 1) = 7o (<fe) + Sf ( 93 ) - 5? (q 3 ) , (B7) 
leads to 

7ife + l) =71(93)- (B8) 

Since only the differences between the Berry phases of 
different states matter, 71 (g 3 ) are defined up to an ar- 
bitrary overall phase. This overall phase can be chosen 
such that, for stage I, 71(93) = 0. 

The Hamiltonian at the end of stage II, H24, can be 
written as 



#24 = - M e^ X 4 t X2T 



h.c 



h.c. 



(B9) 



In the second line, we have used the explicit form of X2,4t 
in terms of the spin and charge operators. Writing the 
Hamiltonian in the basis of eigenstates of e %7T< ^ 2 , we get 

2m -1 

#24 = -M J2 ^ 2i+l " q2 \l2 + l){q2\+h.c. (BIO) 

92=0 



H- 



23 



1*23 I 



E 

92=0 



\92){d2 + A+h.c. 



(B3) 



The above form can be derived from the relation 



ator for e l7r< ^ 2 . The Hamiltonian (BIO I can be thought of 



is a raising oper- 



as an effective tight-binding model on a periodic ring of 
length 2m with complex hopping amplitudes. Note that 
the total effective flux through the ring is given by 



(R23 



2mif2 



(B4) 



In order to diagonalize i?24, we perform a gauge trans- 
formation to a new basis \£ q2 ) defined as 



192, 



16 



q-2 1 



(Bll) 



This transformation is designed such that, in the new 
basis, the phases of the hopping amplitudes are uniform. 
The Hamiltonian takes the form 

2m 

H 2i = - 1*241 E e 4 (^-^)|^2+i)(e 92 l +h.c. (B12) 

92=0 



1G 



which is easily diagonalized in the basis of plane 
waves. One can verify that, defining |£ p ) — 



/T^T l^q 2 =0 C IS©/) 



H 2i \i P 



-2 |t24| COS 



:(p- 



<^24 



14). (B13) 



Therefore, for — fern < </?24 < — (^in + 1) where km is 
an integer, we get that the ground state occurs for p — 
fciiimod (2m). Using Eq. (Bll), one can express the 



ground state in terms of the original basis states: 



|^/( ? 3)) - |$ n («0) 



'2m 



2 2m-l 

£ 

92=0 



(B14) 



Applying the symmetry operator En 
sides of dBlL pit]}, we get 



= e ~™ Sl to both 



Sf(qs) = -—fen, 4 I (g 3 ) = 0. 



(B15) 



Therefore, solving the recursion relation (Eq. B7 1 and 



choosing a gauge such that 711(93 = 0) = 0, we get 



7ii(<?3) 



(B16) 



For the last stage of the evolu tion, Appl ying Eni 
e -inQ 2e i7T§ 3 to both sides of Eqs ^ (B14|,(B1 1, we get 



(93) 



-(sa-fci + i). 

(B17) 



Inserting this into Eq. B7 and solving for 7111(53), we 
obtain 



7111(93) 



2m 



(93 ~h- km) 2 



(B18) 



The total Berry phase 7 = = 71 
unimportant overall phase, is 



7n + 7iii- up to an 



7(93 



2m 



(93 - kf 



(B19) 



Here, k = k\ — k\\ + km. Note that while 7 Q depend on 
our various gauge choices for the basis of the eigenstates 
of H 12, H 23, and H24, the Berry phases of the entire path 
(Eq. B19) does not depend on these gauge choices. 



Appendix C: Topological Protection of the Braid 
Operations 

1. Independence of microscopic details 

In any physical realization, one would not be able to 
control the precise form of the Hamiltonian in each stage 
of the braid process. It is therefore important to discuss 
to what extent the result of the braiding depends on the 
details of the Hamiltonian along the path. Below, we 





Figure 5. Braiding paths in Hamiltonian space, (a) Path 
Pi , for which we compute the braiding adiabatic evolution op- 
erator explicitly, (b) A different path P2, whose Hamiltonians 
at the intermediate stages are assumed to be adiabatically 
connectable to those of Pi. P3 is a path equivalent to P2, 
in which each intermediate Hamiltonian of P2 evolves to the 
corresponding Hamiltonian of Pi and then back. 



argue that the braiding is "topological", in the sense that 
it is independent of these precise details. 

Let us begin by noting that the evolution operator de- 
scribing the full braiding process depends only on: 

1. The initial and final Hamiltonians at each stage; 

2. The symmetries of the Hamiltonian at each stage; 

3. The fact that the ground state degeneracy through- 
out the process is fixed, such that the evolution can 
be considered adiabatic. 

One can see that the precise details of the time- 
dependent Hamiltonian during the braiding process are 
unimportant for our derivation of the evolution opera- 
tor in Sec. I VI I Note that we have never used the exact 
form of the Hamiltonian during the path to determine 
the evolution operator. 

In order to make this argument more formal, let us de- 
fine Pi as the closed path in Hamiltonian space, H12 — > 
H23 — ► H24 — > H\2 , for which we computed the evolution 
operator. (Pi is summarized in Table [I]) Suppose that 
we replace Pi by a different, "realistic" path P 2 , defined 
as H12 — > H23 — > H 2 4 — > H12, which has the same sym- 
metries as those of the original trajectory in each stage 
(Table [I]). Pi and P2 are represented in Fig. [5] a and b, 
respectively. We assume further that the Hamiltonian at 
the end of every stage of P2 is adiabatically connectable 
to that of the original path Pi, e.g. #12 and P^i2 are 
adiabatically connectable, etc. We argue that the adia- 
batic evolution associated with P2 is unitarily equivalent 
to that of Pi. To show this, consider the modified path 
P3 shown in Fig. [5] 



H12 — > H\2 — > H23 ~> H23 — > P^23 

— >P^24 H24 — > H 2 4 — > H\2 Hi2- 



(CI) 



Clearly, Eq. (CI I can be viewed as a deformed version 



of the original trajectory Pi, in which the intermediate 
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Hamiltonian during each stage is deformed relative to 
the original trajectory of Table ^ Since the intermedi- 
ate Hamiltonian in every stage of trajectory (CI I has 



a) 



the same symmetries as those of the original trajectory, 
the analysis outlined in the previous section shows that 
the e volution operator representing the overall trajectory 
(CI I is where to is global phase factor. 

On the other hand, we can consider P3 starting and 
ending with H12, as the "realistic" trajectory P2. Since 
for example, the step H23 —> H23 is "undone" by the next 
step H23 — > H23, Pi is unitarily related to the P3 by 
e lu) VU34V\ where V represents the evolution from H12 
to Hi2- In essence, the matrix V relates the eigenstates of 
the "realistic" initial Hamiltonian H12, to those of H 12 - 
We conclude that the adiabatic evolutions correspond- 
ing to paths Pi and P2 are physically equivalent, and 
therefore the braiding operation is robust to changes in 
the path in Hamiltonian space, as long as the conditions 
1 — 3 listed above are met. 
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Figure 6. Constrictions in an FTI bar. We view the FTI 
as a combination of two FQH droplets with opposite spins 
and opposite filling fractions. The arrows show the direction 
of propagation of the edge states of the two FQH bars. The 
red stars mark the position of the SC/FM interfaces on the 
edges, (a) Constriction created by a gate potential, which 
acts on both spin species. In this case, quasi-particles of both 
spin species can tunnel across the constriction, (b) Applying 
a local Zeeman field splits the spin up and spin down edge 
states, such that only one spin species can tunnel across the 
constriction. 



2. Symmetries of the Hamiltonian during the 
braiding process 

Next, we discuss the symmetry requirements in every 
stage in more physical terms. The topological stability of 
the braiding operation depends crucially on the symme- 
tries of the Hamiltonian throughout the different stages 
of the braiding operation. We now argue that these sym- 
metry properties are largely independent of the micro- 
scopic details of the Hamiltonian in each stage. This 
is since the definition of the braid operation only con- 
tains information regarding which interfaces are brought 
in proximity at each stage. For instance, any Hamil- 
tonian trajectory corresponding to this braid operation 
only couples interfaces 1, 2 and 3 during stage I (see 
Fig. [2]). Any such Hamiltonian necessarily commutes 

with e l7rS:i and e I7r< ^ 3 , independently of its microscopic 
details. For example, adding terms such as higher pow- 
ers of X2,tXi ft retains these symmetries. Likewise, terms 
representing direct coupling between interfaces 1 and 3, 
such as powers of X3,tXi f can be added, as long as they 
are absent at the beginning and end of stage I, when 
interfaces 1 and 3 are far apart. 

A similar statement can be made for stage II: as long 
as interfaces 5, 6, and 1 remain decoupled throughout 
the evolution, the Hamiltonian necessarily maintains the 
same symmetries as those in in Table |T] regardless of the 
microscopic details of the process. 

The symmetry requirement in stage III requires more 
care. At this stage, interfaces 1,2 and 4 are coupled. 
Cruc iall y, we note that the commutation relations in 

= for any j, k 7^ i. There- 
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fore, as long as we allow tunneling of only spin-up parti- 
cles between interfaces 1,2 and 4, we are assured that 
X3.\ commutes with the Hamiltonian. It follows that 
X5 4.X3 1 = e -mQ-2 e inS 3 ^ symmetry operator required 



in stage III, also commutes with the Hamiltonian. Again, 
this symmetry would be maintained independently of 
the exact form of the Hamiltonian, as long as it obeys 
the above restriction. The physical reason behind this 
symmetry is clarified by noting that transferring n up- 
spin quasi-particles to interface 4 (from either 1 or 2) 

changes e"^ 2 e l7r( -® 2+ ™ ) and e in§3 e "^ 3+ ^\ leav- 
ing e - l7r Q2 e "iS 3 j nvar i an t. 

We now see why it is crucial, in order to allow for the 
braiding operation, to have only one species of quasi- 
particles tunnelling between interfaces. If quasi-particles 
of both spins are allowed to tunnel, e _I7r( 5 2 e I7r ' S3 would 
cease to be a good symmetry - in fact, the symmetry 
of the Hamiltonain is lowered, and the ground state de- 
generacy is reduced from (2m) to just 2, violating the 
adiabatically of the braiding process. This is a special 
property of the fractional (m > 1) case; for m = 1, there 
is no difference between up and down quasiparticle tun- 
neling. 

We note that the restriction to single species tunneling 
was unnecessary in stages I and II, which retain the same 
symmetries even when both spin species are allowed to 
tunnel. Moreover, if we allowed only spin down quasi- 
particles to tunnel between 2 and 4, there would be an 
alternative symmetry operator e %7T< ^ 2 e l * S;> at stage III. 
The braid operation with this type of coupling would 
yield a unitary operator of the same form found in Sec |VI| 

The restriction of single spin species tunnelling can be 
met in the different realizations of the model analyzed 
above. First, consider the realization using a fractional 
quantum Hall liquid, in which an insulating trench sep- 
arates two counter propagating edge states (Fig. [TJd). 
In this realization, the labels spin up and down indi- 
cate whether the quasi-particles are on the inner or outer 
edge respectively. When we deform the system in order 
to put interfaces in proximity, we must specify whether 
this deformation shrinks the inner or outer droplet of 
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the quantum Hall liquid. Suppose we shrink the inner 
droplet. Then quasi-particle tunnelling between inter- 
faces proceeds through the inner droplet. Therefore in 
this case, quasi-particles can only tunnel from the inner 
edge at one interface location to the inner edge at an- 
other. Only electron tunneling processes are allowed be- 
tween the interfaces on the outer edge. However, electron 
tunnelling from the outer and inner edge are equivalent, 
as these are related by cooper-pair tunnelling or "spin 
flip" operators. To conclude, in this realization, choos- 
ing whether to deform the inner or outer quantum Hall 
droplet selects which spin species of quasi-particles are 
allowed to tunnel between interfaces. 

Let us now consider the realization of the system on the 
edge of a Fractional Topological Insulator. Suppose that 
one can apply either ordinary gate potentials, or Zeeman 
fields in the z direction (by coupling to a nearby ferro- 
magnet polarized along z). which act as opposite gate 
potentials for the two spin species. Then there are two 
ways of coupling two interfaces, depicted in Fig. [6] One 
can either create a constriction in both spin species by 
applying an appropriate gate voltage (Fig. [6^), which 
allows quasi-particle tunnelling of both spin species be- 
tween the interfaces across the constriction, or create a 
constriction for one spin species only, e.g. spin up (Fig. 
(6j)); in which case only that spin species tunnels. Note 
that in the latter case, we have split the spin up and 
spin down edge states into two counter-propagating edge 
modes, which become gapless. However, if the length 
of the split region is L, there still is a finite-size gap of 
the order of vp/L, where vp is the Fermi velocity on 
the edge. The tunnelling of quasi-particles of spin up 
across the constriction, on the other hand, is enhanced 
by a factor of the order of exp (L/£) relative to that of 
spin down, where £ is the correlation length in the bulk. 
Therefore, the tunnelling of spin up quasi-particles can, 
in principle, be enhanced parametrically without reduc- 
ing the gap considerably. 



U13 = e % — 



U 23 = e 1 * 




e-^ k +''f e " rfes . (Dl) 



fe=0 



Here, we have used the expansion of the braiding matri- 
ces in terms of the spin and charge operators and their 
harmonics. 

The Yang-Baxter equations state that 



U12U23U12 — U23U12U23. 



(D2) 



This relation can be understood pictorially, as shown in 
Fig. [3] Inserting Eqs. (Dll into the left hand side of 
(D2|, and using e ™Q e ™& = e -i^ e iwS e inQ^ we get 



u 12 u 23 u 12 = ]T 

kiMM (2m) 



ki,k 2 ,k 3 



(2m)- 



e i7rfci Q e iirk 2 S gi7rfc 3 Q 



_ giTtk 2 S giir ( fci +fc 3 ) C 



(D3) 



The sums over ^1,2,3 run from to 2m — 1. Changing 
variables k% — > k\ + &2 , and k 3 — > k 3 — k\ , 



, p -i^[k 2 2 + (k?-ki) 2 ]+i^ 

u 12 U2 3 u 12 = J2 1 3 

ktMM ( 2m ) 2 

e — '5m (ki+kl+kl^+i '^- 

kuk 2 ,k 3 (2m) 5 



gi7rfc 2 Sgi-n-fc3 Q gi-7rfci S 



1,K 2 ,K 3 

e i^s 2 e i^Q* e i^s- = U23Ul2 u 23 , (D4) 



Appendix D: Yang-Baxter equations 



Here, we verify that the unitaries representing braid- 
ing of two neighboring interfaces by tunneling of quasi- 
particles satisfy the Yang-Baxter equations. Imagine 
that we start from three consecutive interfaces, 1, 2, and 
3, shown in Fig. [3] The segment between 1 and 2 is a 
superconducting (SC) segment, and the segment between 
2 and 3 is a ferromagnetic (FM) segment, e 171 ® and e l7rS 
are the charge and spin operators acting on the SC and 
FM segments, respectively. In terms of these operators, 
one can express the unitary matrices that correspond to 
braiding (1,2) and (2,3): 



In the second line we have commuted e l7rfclS with 
e iirk 3 Q _ This establishes the Yang-Baxter equation for 
the braiding matrices (Dll. Conjugating Eq. (D2| by 



e nrQ L k le ~nrS R k 2 ^ wnere ^Ql correspond to the charge in 
the SC segment on the left of e i7r< 3 ) anc l e ^S R ^ Q ^ ne spm 
in the FM segment on the right of e l7rS , we get Eq. (27 1. 



Appendix E: More on the Braiding and Topological 
Spin of Boundary anyons 



In Sec |VII| we have derived the representation of the 
braid group using an analogy to braiding properties of 
anyons in two dimensions. The derivation proceeded us- 
ing two important assumptions, and below we explain 
why these assumptions actually follow from properties of 
two dimensional anyons. 
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(i) 



(ii) 





Figure 7. Properties of particle exchanges (i) Two con- 
secutive exchanges, first between particles a and b, and then 
between particles b and c. (ii) Particle a is exchanged with a 
particle d, which is the particle resulting from fusing b and c. 
For two dimensional anyons, a proper definition of the braid 
matrix between any two particle types establishes an equality 
between (i) and (ii). For the case a = b — c = X and d = q, 
we only use the property that the two braids in (i) move the 
label q one segment to the left, as shown in (ii). For the case 
where a = qi,b = qi,c = qz and d = qi + qz mod 2m, we 
show that the two sides are indeed equal 



1. Properties of the particle exchanges 

Consider a two dimensional theory, in which particle of 
type a is exchanged first with particle of type b and then 
with particle of type c. The operation should depend 
only on the type of particle a, and the total topological 
charge of particles b and c. The operation should not 
be able to distinguish the finer splitting of the combined 
charge into the charges b and c. This condition can be 
summarized pictorially in Fig [7J 

Importantly, the conditions summarized in Fig [7J hold 
also in our one dimensional system, for the representation 
derived in Sec. |VII| as we shall show below. The impor- 
tance of these conditions to our derivation in Sec. I VIII is 
threefold, since it can be used to (i) show that the pro- 
cess in Fig. |4|c) indeed exchanges the charges of the two 
segments; (H) show that the TS of the q±, g 2 composite 
is equal to the TS of q\ + q 2 mod 2m; (iii) show that the 
Yang-Baxter equations hold for a braid of any two parti- 
cle types. Therefore, the only assumption necessary for 
the derivation presented in Sec. |VII| is that the condition 
summarized in Fig. [JJ holds. 

In the following, we shall study three important cases, 
for which we verify explicitly that the conditions in Fig.[7J 
hold in our one dimensional system. These three cases 
are used to verify the assumptions (i) and (ii) above. 



a. a — b = c = X,d — q 

For concreteness, Consider four neighboring inter- 
faces X1-X4, flanking alternating FM, SC and FM seg- 
ments. The operators corresponding to the charges in 

= e^Q, X 2tX t i T = 



the different segments are Xst^i']' 



We shall also use xv\X 



Consider an initial state of the system 
which is an eigenstate of e m ® with eigenvalue e l ™ q , and 
of e Z7rS i+ s 2 w ith eigenvalue e l ™ Stot , 



k,stot) 



2m-l 



- qn 



1 2m n 

n— 



Si 



-n, s 2 = stot + n) 



(El) 

where s\ and S2 correspond to e I7rSl and e l7lS2 . Note that 
such a state is also an eigenstate of X4TXi^- 

Following Fig. [7] consider two braid operations first be- 
tween X\ and X2 and then between X2 and X3. Using 
e l7r( 2|si, s 2 ) = \si — l,s 2 + 1), and the Fourier represen- 
tation of the braid operators, the resulting state is 



(2m) e 

2m- 1 

E ■ 

n,p,fc— 



^[-(p 2 +fe 2 )+2< 2 n-2np]| _ 



(E2) 

k, stot + n + k) 



Denoting I — n + k we arrive at 

(2m) e - l 5;y 12 (Q)[/ 2 3(^i)k,Stot) 

2m -1 

= e i 2^h((P+ n ) 2 +' 2 )+ 2 («+') n ]| - I. 

71,p,i— 

and therefore 



(E3) 



Stot 



•0 



(2m) e-^U 12 (Q)U 23 (S 1 )\q ) s tot ) 



(E4) 




Si 



q, s 2 = stot - q) 



We see therefore, that the two consecutive exchanges are 
equivalent to moving the charge q one segment to the 



left, i.e., e l »» becomes the eigenvalue of X2tXit = e 

That is exactly what is indicated in Fig. |jj (b). The 
state is also multiplied by a phase which depends on the 
gauge choices for the different basis^. Importantly, note 
that we could have chosen to fuse X\ with a different 
interface Xi (as long as it is not between 1 and 3), to 
form a charge XifXif) which would again commute with 
X3t%2f identical analysis to the above would yield 
the same result, with Xi-f replacing Xi-]-- 



b. Exchange of two q 's 

Using the above, we would now like to show that 
the four exchanges depicted in Fig. [4] (c) indeed corre- 
spond to exchanging two charges q\ and q2- Consider 
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four neighboring interfaces Xq-X 3 , flanking alternating 
SC, FM, SC segments, where the initial state is \qx,q 2 ), 
corresponding to an eigenstate of XitXo^ = e l7rQl and 

Indeed, using the above results, we see that performing 
the exchanges X\ with X2, and then X2 and X3, would 
result in an eigenstate of X2fXi-t- = e ™ S with eigenvalue 
e 1 ™ 92 . Now, performing the exchanges Xq with X\, and 
then X\ and X2, and again using the previous results, we 
see that the resulting state is an eigenstate of XitXof — 
e l7V Qi with eigenvalue e 1 ™ 92 . Therefore, it is an eigenstate 
of X3-|-X2f = e l7r ® 2 with eigenvalue e ! ™ ?1 (the total charge 
in the two segments is preserved). Therefore, the two 
charges have been exchanged by the sequence depicted 
in Fig. |4](c). 

We shall now explicitly calculate the phase factor re- 
sulting from an exchange of two charges. Let us denote 
the operation at point by Uq i q 2 . We would now like to 
verify that 



U^^\q u q 2 ) = U 23 (S)U 34 (Q 2 )U 12 (Q 1 )U 23 (S)\q 1 ,q 2 ) 

(E5) 



= e^«^\q 2 , qi ) 



and find the abelian phase 4>{qi,q 2 ) associated with this 
exchange. The sequence of the four braid operation in 
the above equation corresponds to Fig [3] (c). Using 

s z ' KS \qi,q2) = \qi — 1 7 <Z2 + 1) and the Fourier representa- 
tion of the braid operators, we arrive at 



e l ^( n2+p2 ) +{qi - n '> 2 + iq2+n) ^\q 1 ~n~ P ,q 2 + n+p) 

n,p—0 

(E6) 

Denoting I = p + n we arrive at 

2m-l 

e l ^- l2+2nl ) +2{q2 - qi)n+q2+q2 }\ qi -l,q 2 + l) (E7) 



,1=0 



The sum over n forces I 



qi — q 2 . which therefore gives 

(E8) 



e^ q ^\q 2 , qi 



Therefore, Eq. (E5l hold, and the abelian phase for in- 



terchanging two charges qi and q 2 is just e 1 ™ 9192 , which 
is consistent with the topological spin for the charges, 



Eq. (291 



first have an exchange of the charge of the first and 
second segments, resulting, according to the above dis- 
cussion, in the state e 1 ™ 9192 \q 2 , qi, (J3), followed by an 
exchange of the second and third segment, resulting in 
e l ^^ qiq2+qiq ^\q 2 , q 3 , q^ . Clearly, that is exactly the re- 
sult of the operation on the left side of Fig. [7] in which 
the total charge q 2 + q 3 is exchanged with qi and then 
split again into the two segments. 



2. Consistency check 

In this section we would like to verify that the topo- 
logical spin of the interfaces is indeed well defined, i.e., 
the phase acquired by the operation defining the topolog- 
ical spin is independent of the initial state of the system. 
Consider the procedure in Fig. [4] (e) . Interface 1 flanks 
the left side of a SC segment, which corresponds to e l7r( 3. 
The initial state can be taken as an eigenstate of e nCl 
with eigenvalue q. A FM segment is nucleated to the 
right of it. The FM segment is nucleated with spin zero. 
Therefore, the state after the nucleation of this segment 
is 



1 



2m-l 



'2m 



X] \qi = q-n,q 2 = q + n) (E9) 



n=0 



|#7,s = 0)) = 



As explained in Sec. |VII[ the braid operation between 
interfaces 1 and 2 does not change the charge of the SC 
segment between interfaces 1 and 2, the charge between 
3 and 4, and therefore the total charge q (numbering 
increases to the right). Therefore, 



2m- 1 



cr 12 hK 9 ,«)) = -L 53 

V2m r, 



e^ q -^\ qi = q -n, q2 = q + n) 



n=0 



(E10) 

where we keep the phase function ip(q) completely gen- 
eral. We now project on the subspace with e l7rS = 1 in 
the FM segment, i.e. we apply the projector 

n 5=0 = 53ivfes = o))(^( s = o, g )| (eh) 



were the identity operation is implicitly assumed to act 
on all other degrees of freedom. Applying the projection 
yields (up to normalization), 



c. a = qi, b — q2, c — qs d — q2 + q L i mod 2m 

Consider three consecutive SC segments, and an ini- 
tial state \qi,q 2 ,q 3 ), corresponding to an eigenstate of 
e"^ 1 ,e^ 2 ,and e 1 ^ 3 . On the left side of Fig. we 



2m- 1 



? X = eMq ~ 



(E12) 



n=0 



Importantly, Ox does not depend on q, as the sum runs 
over all possible values for charges. 
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